Model study of dissipation in quantum phase transitions 
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We consider a prototypical system of an infinite range transverse field Ising model coupled to a 
bosonic bath. By integrating out the bosonic degrees, an effective anisotropic Heisenberg model is 
obtained for the spin system. The phase diagram of the latter is calculated as a function of coupling 
to the heat bath and the transverse magnetic field. Collective excitations at low temeratures are 
assessed within a spin- wave like analysis that exhibits a vanishing energy gap at the quantum critical 
point. We also consider another limit where the system reduces to a generalized spin-boson model 
of two interacting spins. By increasing the coupling strength with the heat bath, the two-spin 
wavefunction changes from an entangled state to a factorized state of two spins which are aligned 
along the transverse field. We also discuss the possible realization and application of the model to 
different physical systems. 
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I. INTRODUCTION 

The twin (and apparently disjoint) topics of quantum 
phase transition (QPT) and quantum dissipation (QD) 
have seen a great upsurge of activity in recent years. The 
QPT is of significance in many areas of contemporary in- 
terest in the Condensed Matter, such as metal-insulator 
transition, quantum magnetism, ferroelectricity, super- 
conductivity and in the general question of coherence 
and quantum computation^. On the other hand, QD is 
ubiquitously present because of environmental influences 
on otherwise unitary evolution of quantum systems-. 
Dissipation, though considered a pest, needs to be under- 
stood (and tamed) , in order to tackle decoherence effects 
in quantum many body systems. Our aim in this paper is 
to analyze the combined presence of these two seemingly 
disparate phenomena of QPT and QD via simple model 
systems. Our hope is to elucidate on the irreversible ef- 
fects near a quantum critical point(QCP) because of dis- 
sipative interactions with the environment. 

The simplest model of a QPT is an Ising system of 
coupled spins (which are taken to be polarized in one di- 
rection only, say the z-axis), subjected to an additional 
('magnetic') field T along, say the x-axis. The latter cou- 
ples to the x-component of the spins which, because of its 
non-commutativity with the z-component, triggers quan- 
tum dynamics in the system. The resultant 'Transverse 
Ising Model (TIM)' is a prototype for analyzing the con- 
flicting presence of 'order' along z-axis induced by the 
Ising coupling and 'disorder', caused by the tilting of the 
spins to the x-axis by T. The net result is the occurrence 
of a QCP at temperature T = in the V — T phase di- 
agram as r is increased to a critical value r c , when the 
system transits from a ferromagnetic to a paramagnetic 
phase (see Fig.l)^. We will be interested in analyzing 
questions such as what is the analogue of 'phase ordering' 
in classical systems (see^) when the quantum mechani- 
cal system of Fig.l is subjected to a sudden 'quenching' 
along the T-axis across r c , maintaining the temperature 
at T = 0. Because a change in the value of T (from above 
T c to below) will inevitably lead to an irreversible transi- 



tion from one equilibrium configuration to another, dis- 
sipation needs to be dovetailed to the discussion. Again, 
a straightforward and widely studied model of quantum 
dissipation, popularized in recent years by Ford et al^, 
and Caldeira and Leggett£, is the one in which the quan- 
tum system at hand is linearly coupled to a bosonic bath. 
TIM is an effective model which has application from 
solid state materials of rare earth magnets to a coupled 
Josephson arrays^, where dissipation arises in a natural 
way. Our model Hamiltonian can then be written, 



^=J2 s t9k{bl + b k ) + J2^kb{b k , (1) 



where iV is the total number of spins and gk is a coupling 
constant. When gk = and the Ising coupling is 
treated in the mean field approximation(MFA), we obtain 
the phase diagram, depicted in Fig.l. 

Several limiting cases of Eq(p]) have received attention 
in recent literature. If Jij — and the term sf in the cou- 
pling with the heat bath is replaced by sf, Eq([T]) yields 
the celebrated spin-boson model of quantum dissipation 8 . 
When the Ising interaction is replaced by a Zeeman cou- 
pling with an external field ft on a single (N — 1) spin, 
the model in Eq([T]) describes low-temperature dissipa- 
tive quantum tunneling in an asymmetric double well^— . 
Additionally, if h — 0, one has a symmetric double- well, 
tunneling in which can be impeded, leading to localiza- 
tion, when the coupling with the bosonic bath exceeds a 
certain critical valu o 12 ' 13 . This case is also relevant for 
a Kondo impurity of spin one-half (described by s), in 
interaction with a conduction electron bath, which can 
be modeled in terms of bosons as far as the electron-hole 
excitations near the Fermi surface are concerned^. For 
r — 0, this model can be viewed as an generalization of 
the Dickc model of 'superradiance' where the Ising term 
includes an additional atom-atom interaction^. 

One other important application of Eq([T]) ensues in the 
case wherein the range of the Ising interaction is infinite, 



2 



- 








FM \ 


PM 



FIG. 1: Mean field phase diagram of a TIM in which the 
abscissa represents the transverse field F whereas the ordinate 
depicts the temperature T. When r = 0, we have the usual 
Curie transition at T c . On the other hand at T — 0, there is a 
QPT at r c from a ferromagnetic(FM) to a paramagnetic(PM) 
phase. 
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a constant. This situation is the one in which the MFA 
to the Ising model (in the absence of the transverse field 
r) becomes exact and will occupy much of our attention 
below. Writing the total spin components as 



E' 



5. 1! 



(3) 



we have from Eq([T|) and Eq(J2J) , 
J 



H 



2N 



k v k 

(4) 

If we leave aside the coupling term gi~, the spin- 
Hamiltonian in Eq((4]) represents a 'Molecular Magnet' 
characterized by a single-ion anisotropy energy J/N, and 
subjected further to a transverse field T—. This prob- 
lem has been investigated in great detail in recent years 
in the context of 'macroscopic magnetization tunneling' 
when the value of S can be pretty large such as S = 10. 
The additional coupling to the bosonic bath when g^ is 
switched on, enables us to treat the effect of dissipation 
on this tunneling behavior—. The preceding remarks 
then underscore the versatility and relevance of the model 
Hamiltonian H in Eq([T|) in a variety of applications to 
current topics in the condensed matter. In the sequel 
we shall analyze diverse mean-field and low-temerature 
properties of H, keeping the underlying QPT in mind. 
With this background to the formalisms developed here, 
the paper is organized as follows. 

In section II we first carry out a unitary transforma- 
tion on the Hamiltonian in Eq([T|), that has been bor- 
rowed from the literature on polaron-physic o i ' 18 . This 



transformation enables the original coupling constant 
5^ (proportional to a variational parameter /j.) to be el- 
evated to an exponential function, thus facilitating an 
analysis that works even in the regime of strong coupling 
to the bath. The transformed Hamiltonian, rewritten in 
terms of fk, is then treated in section III in the MFA 
to evaluate the associated density matrix from which the 
Helmholtz free-energy and an equation of state for the 
magnetization can be calculated. In section IV we fo- 
cus our attention to spin wave like excitations near the 
absolute zero of temerature. The section V is devoted 
to a novel aspect of entanglement, important in the con- 
temporary issue of quantum information process, in the 
context of two coupled Ising spins. Finally in section VI 
we present a brief summary. 



II. EFFECTIVE SPIN-HAMILTONIAN 

In this section we derive an effective Hamiltonian of the 
spin system, starting from Eq(p}. For this, we integrate 
out the bosonic degrees of freedom of the heat bath. Here 
we adopt the usual definition of the 'effective partition 
function' Z e ft of the spin system^, 



J eff 



-piHs+HB+Hr 



/Tr B e 



-BHz 



(5) 



where H s is the Hamiltonian of the 'transverse field Ising 
model', Hb describes the non- interacting bosonic degrees 
of the heat bath and Hi is the interaction term with 
the heat bath. To decouple the spin-system from the 
bosonic heat bath, we subject the original Hamiltonian 
to a unitary transformation, 



H = UHU~ 



(6) 



where, 



U = exp 



(7) 



At this stage we treat as a variational parameter 
which can be determined from the minimization of the 
free-energy of the total system. As a special case, if we 
consider a non- interacting spin system by setting Jy = 0, 
we notice that the total Hamiltonian can be diagonalized 
by the unitary transformation given in Eq.© with fk = 
gk/hujk- Motivated by this observation, the above men- 
tioned variational method has been applied to a single- 
spin(V = 1) 'spin-boson' model, which successfully cap- 
tures the 'Kondo' like localization transitio n 18 ' 20 . 

Following the unitary transformation in Eq.(J7]) on 
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Eq.©, 
H = 



2N l 



N 

T 



+ S 2 Z + -(S 2 Z - S 2 y )(cosH-2F) - 1) 



--(S z S y + S y S z ) sinh(-2F)] - TS X 
^Z(bk + b{ ) (g k - huj k f k ) 



S 



(8) 



III. PHASE TRANSITION OF QUANTUM 
ISING MODEL: N -> oo 



In this section we consider the case wherein the number 
of spins N — > oo, which corresponds to an 'infinite range' 
quantum ising model in the thermodynamic limit. For 
the corresponding classical model it is known that the 
MFA is exact in this limit. 



A. Mean-Field approximation 



where, F — J2 k ^Tiy^fc ~~ ^fe)> and S x>z are components 
of the total spin as given in Eq.Q. Now we approxi- 
mate the total density-matrix as a direct product form 
P = Ps&) pb, where p s and p B denote the density matri- 
ces of the spin system and free bosons respectively. After 
integrating out the bosonic modes of the system, we ob- 
tain the effective Hamiltonian H e f / of the spin system: 



Tr B H = H, 



eff 



(9) 



where Fb is the free energy of noninteracting bosons and, 
H eff =-K z S 2 z -K y S 2 y -K x S 2 x ~rS x . (10) 



The Hamiltonian above describes a fully anisotropic 
Heisenberg-model in the presence of a magnetic field 
along the x-axis, when it is written using Eq. ((3]) . Ef- 
fective coupling strengths are given by, 



K x 



K, 



k 



_J_ 

AN 
J 

4JV 



1 - exp{-—J2fkC0th(hu k /2k B T)} 



l + exp{-—J2f k 2 coth(tu» k /2k B T)} 



1.1) 



Going back to the notation of total spins (see EqQ) we 
see that the terms proportional to K x and K y generate 
new physics in the context of molecular magnets. From 
minimization of the free-energy of the total system with 
respect to the variational parameter f k , we obtain, 



fk — 



l9k{S 2 x )/N 



[2^^ + j£{{S 2 ) - (S 2 )) coth(huj k /2k B T)} ' 



where 



K = exp 



N ^ 

k 



/ fe 2 coth(^ fc /2fc B r) 



(12) 



(13) 



Thermodynamics of the spin-system described by the 
effective Hamiltonian Eq. ([T0|) can be obtained by solv- 
ing the self-consitent equations described in Eq (TT2")) and 
Eq (fT3"|) . as described below. 



Within the MFA the phase diagram of the above model 
can be sloved analytically. Here we assume that the den- 
sity matrix of the spin system can be written as product 
of the density matrix of each spins p s — Yii&> Pi- The 
density matrix of each spin can further be expressed as, 

1 



Pi 



;(i 



(14) 



where the order parameters are (sf ) = m/2 and (sf ) = 
x/2. From the selfconsistency equation for f k Eq (Ti"2"|) we 
obtain, 

Jm 2 K 

■coth(hw k /2k B T)+ X 2 f^k -(15) 



fk 



X 2 9k 



2N 



0, in the limit of N — > oo, 



If we assume J2k 9kl^' u - 
we have K — > 1. 

The free energy of the system is given by, 



F 



JNm 2 



y N\ Nx 2 



s 



k 



(2/fe3fe ~ fr^kfk) 



-k B TN 
-F B 



-(1 + ln(l + + -(1 - ln(l - - In 2 



(16) 



X 2 and F B is the free energy of the 



where £ = \fri 
bosons. Minimizing the free energy per particle F/N 
with respect to the parameters m and x we obtain: 



r 

' 2 



Jm 
2 



kTm . 
kTx 



In 



2£ 



In 







0, 



(17) 
(18) 



where A = J^kdk/^k- The ferromagnetic phase is de- 
fined by to ^ and the phase boundary can be obtained 
from 

jr 



T/(J/2 — A) = tanh 



4fc B T( J/2 - A) 



(19) 



A three-dimensional plot of the phase diagram is shown 
in Fig. 2. As expected of course, for A = 0, Eq. (Tig)) 
yields the usual equation of state and the phase diagram 
depicted earlier in Fig. 1. The latter also bears out the 
expectation that as the strength of the coupling to the 
heat bath A increases, the QCP on the T-axis is sup- 
pressed. 
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FIG. 2: Phase boundary of ferromagnetic phase of a TIM 
coupled to a heat bath (see Eq. (Tl9)) '). Temperature T, trans- 
verse field r and strength of dissipation A are measured in 
units of J. The region under the curve represents the ferro- 
magnetic(FM) phase. 



B. Semiclassical analysis 

In this subsection we calculate the effective partition 
function of the spin-system semi-classically and establish 
the validity of the variational method and MFA for a 
large spin system (N — > oo). We can decompose a direct 
product of N spin 1 /2 systems into a direct sum of total 
spin S, where the magnitude of the spin varies from N/2 
to or 1/2 (depending on even or odd iV). In terms of 
the total spin operator S 1 , the partition function can be 
written as 21 



Z = ^2 Tr B£ 



N/2 

^ = ]T P{S)Tr s Tr B 

5=0,1/2 



(20) 



where P(S) is the number of ways the total spin S can 
be formed 2 ^ and is given by 



P(S) 



(2S+1)N\ 



(N/2- S)l(N/2 + S+l)V 



(21) 



For large S we can calculate the partition function of 
spin systems classically 2 ^ and replace the trace by the 
integral, 



Tr, 



2S+1 

47T 



sin 



(22) 



where the components of the total spin S are given by 
(5 sin cos 0, S s'm 9 cos cj>, Scosd). Using the coherent- 
states for bosons we integrate out the bosonic degrees of 
freedom and obtain the effective partition function of the 



spin system as, 



N/2 



z eff = ]T p(s) { ^±^ , sn i 



4ir 

S=0,l/2 

-P[-ih s2 cos2 0~rs sine cos 4>~±s' 



(23) 



The minimum value of the energy is obtained for = 0. 
For large S we introduce two variables r = and x = 
cos 8. In terms of these two variables we can write the 
effective partition function as: 



Z, 



N r / drdxe-P N i'^ x \ 



(24) 



where /(r, x) is the free energy per particle and N r is the 
a constant. Here we have used Stirling's approximation, 

N 

InP(S) « [(1 + r)ln{\ + r) + (1 - r)ln{l - r) - 2ln2] , 

(25) 

which is the entropy in MFA. The free energy per particle 
is given by: 



f(r,x 



J 2 2 Tr rr 

— r x v 1 

8 2 v 



A 



,2 _ - r 2n_ x ^ + 
4 



[(1 + r)ln(l + r) + (1 - r)ln{\ - r)} 



(26) 



Within a saddle point approximation we minimize 
f(r,x) and obtain the following set of equations: 



df(r,x) 
dx 



0. 



ry 1 — x 
df(r, x) 



J/2- A' 



(27) 



dr 



= 0, 



<? - r^/i -x 2 - Ar(l - x 2 ) 



-k B T\n\ 



■ 1+7- 

l-r J 



0. 



(28) 



Defining the magnetisation along z-axis as m/2 = rx/2, 
and using Eq (|2"5)) at the saddle point of the free energy 
function /(r, x) we obtain, 



1/2 



= tanh ■ 



J 



4fc B T 



(J/2 - A) 2 



which checks with the earlier Eg. (119 



r 2 



(J/2-A) 2 
(29) 



IV. COLLECTIVE EXCITATIONS NEAR T = 

In this section we analyze the low temperature spin 
wave like excitations by calculating the low-lying ener- 
gies of the spin system coupled to a heat bath. In the 



5 



previous section we used the total spin representation 
S a = J2i s iai where a = x, y, z, to evaluate the partition 
function semiclassically. Quantum effects become rele- 
vant near zero temperature and consequently, the total 
spin takes a large value S ~ N/2, in order to minimize 
the free energy. In the limit of large S, we can derive 
an effective Hamiltonian describing the quantum fluctu- 
ations at low temperatures. In the Holstein-Primakoff 
representation, the components of total spin S are given 
by, 



S- = bW2S-tfb, 



S+ = ^2S - Wbb. 



(30) 

(31) 
(32) 



where bQr) satisfy bosonic commutation rules. We as- 
sume that the total spin vector S makes an angle 9 with 
the z-axis and its projection on x-y plane makes an angle 
4> with the x-axis. Now we can write: 

S' z = (S - 6 f 6) cos 9 - y/S/2{tf + b) sin 9 (33) 
S' x = v /^/2[(6 t + b) cos 9 cos <j> + i(b-b f ) sin./)] 

+(S- tfb)sin9 cos <f> (34) 

Substituting Eq (l3"4")l in the Hamiltonian Eq. (HJ) , we ob- 
tain, 

H = (-^[(^-^) 2 cos 2 ^+f(& + 6t) 2 sin 2 

- 2W^(£-&tb)(&t + 6)cos0sin0] 
~S 

— {(6 f +b) cos6»cos<?!)-|-i(6- fe^sin^} 



+ (S - b^b) sin 9 cos <f>] x [T - -^=^g k {b 
+ Y,^kbib k + 0{1/S). 



h)] 

(35) 



ij 



After performing a transformation c\. = 
—7==- — q k sin cos <i, we obtain from the above Hamil- 

tonian the classical energy of large spin S (part of the 
Hamiltonian which is proportional to N for large 5), 



H 



JS 2 
2N 



cos 2 9 — TS sin 9 cos < 



N' 



S 2 sin 2 9 cos 2 . 



(36) 

Again, a minimum of Hq can be found for = and 
cos# = for the paramagnetic state. In the ferro- 
magnetic state of spin 5, the minimum energy can be 
achieved for <b = and, 



W- 

N V 2 



A) sin 9 



(37) 



The above equation determing the minimum classical en- 
ergy is equivalent to the equation derived from the sad- 
dle point approximation to the free energy described in 



sec. HB (see Eq (|28| ) . With this choice of 9 and cf>, we 
notice that the terms linear in fluctuation operators (6 
and Cfc) in the Hamiltonian (Eq([3"5j)) vanish. 

Finally the Hamiltonian describing the fluctuations 
can be exactly mapped on to the Caldeira-Leggett 
model 6 describing an oscillator coupled to a heat bath, 



H 2 = etfb + A(6 f + bf 



hu k c\c k (38) 



where different parameters in the above Hamiltonian can 
be written in terms of the spin S and its orientation 9, 



JS 



2XS 



e = — cos 2 9 + —sin 2 6 + T sin 9, (39) 
A = - W sin 2 9 } (40) 

(41) 



!>''■ = 9k\ TV ri,> ' 



Following the work of Ambegaokar and Hakim 2 ^, we can 
diagonalize the above Hamiltonian quadratic in bosonic 
operators by means of the canonical transformations, 

c a = ^ A a pcp + BgpCg, 



4 = ^B a pCp + A a pc\, 



(42) 

where we denote the original bosonic operators as {c Q } = 
{&, Cfc}. The Hamiltonian can be written in diagonal form 
H 2 = y^,a EpCgCf) in terms of a new set of bosonic oper- 
ators. The set of energies Ep describes the low- lying ex- 
citation energies of the many-body system. Equation of 
motion of the operators can be obtained from the Hamil- 
tonian, yielding 



db 
i— 
dt 

A 

*~JT = ^kCk + 5fc(fo f + b), 
dt 

dep . 



eb + 2A(tf +b) + Y,9k(c{ + c k ), (43) 
k 

(44) 
(45) 



Substituting Eq (|42"]) in Eq (|45]) . the matrix elelements of 
the canonical transformation are obtained as, 



EpA p 

-EpBop 

EpA k p 
-EpBf.p 



eA a p + 2A{A p + B p) + J29k(A kp + B k/3 ), 

k 

eB p + 2A(A Q p + B p) + ^g k (A k p + B k p), 



t^kA k p + g k (A k p + B k p), 
fwkB k p + g k (A k p + B k p). 



(46) 



6 



Finally, we obtain the equation determining the exci- 
tation energies Ep, 



4 



4eA 



gpko h 



(47) 



At T = the total spin is S = N/2, and we can define 
two phases according to the orientation of the large spin. 
When r > J/2 — A, all spins are aligned along the x-axis, 
hence cos 9 — solution minimizes the classical energy of 
the system. The QCP is given by, 



T c = J/2 - A. 



(48) 



In the paramagnetic phase the fluctuation operators of 
spin (&T ,b) decouple from the modes of the bath in the 
leading order, since 



9k = 9k cos 9/2 = 



(49) 



and they are coupled in the higher order terms of the 
Hamiltonian which are supressed by a factor of 1/y/N. 
Quantum fluctuations of the paramagnet can then be de- 
scribed by an effective harmonic oscillator 



H 2 = fujpclca, 



with the excitation energy: 



hLo p = v / (A + r)(A + r- J/2). 



(50) 



(51) 



Here we note that although the fluctuations of spins are 
decoupled from the bath modes, the frequency ui p de- 
pends on dissipation. At the QCP, the excitation energy 
vanishes as (r— Tc) 1 / 2 , which is in accordance with mean- 
field behavior. 



V. GENERALIZED SPIN-BOSON MODEL FOR 
TWO INTERACTING SPINS 

In this section we address the issue of quantum infor- 
mation and Schrddinger cat like state by considering our 
model-Hamiltonian Eq.([T]) with N = 2, which represents 
two interacting spins attached to a heat bath. This gener- 
alizes the usual spin-boson model^ by including spin-spin 
interaction terms. 

Following the procedure of unitary transformation and 
then tracing out the bosonic degrees of freedom (as men- 
tioned in section II), we obtain the following effective 
Hamiltonian for the total spin, 



H = -K Z S 2 Z - K v S 2 y - K X S 2 X - TS : 



(52) 



which follows directly from Eq. (1101) . Since at zero temer- 
ature, the triplet state plays an important role, we ob- 
tain the following eigenvalues and eigenfunctions for the 
triplet state: 



e± = ~{K Z + K y )/2 -K x ± J(K Z - K y ) 2 /A + T 2 (53) 




FIG. 3: Variation of the parameter cos 29 of the wavefunction 
|— ) (which is also the 'concurrence' of the wavefunction) with 
the change of the coupling a with the heat bath, for Y/J — 



10" 



f2 c /J = 10 2 (solid line), Y/J 



10" 



n c /j 



(dashed line) and Y/J = 10 , fi c /J = 1 (dot-dashed line). 



±>=^(|tt) + IU» 



sin v± 



V2 

with tan#-t = — (e± 



V2 

K z +K x )/T, and 



(ltt) + Ut», (54) 



£2 = - Ky, 

|2) = ^(|tt)-|U». 



(55) 
(56) 



The singlet state has zero energy. For a system without 
dissipation( when K y = 0, K x = 0) the ground state 
wave function is I— ), with 



cos(26>) 



K z /2 



(57) 



Unlike the Ising system in the thermodynamic limit (for 
N — > oo), this state does not have a net magnetization 
and (S z ) — 0. The magnetization along the x-axis in- 
creases smoothly with increasing the magnetic field T 
and finally for very large T, the state |— ) is factorized 
in two spin states directed along x axis. 

Now switching on the coupling with the heat bath, we 
obtain, 



cos(26>) 



JK/8 



(JK/8) 2 + T 2 



(58) 



For the state R,^ 2 ) = 1,(5 2 ) = cos 2 9,{S 2 y ) = sin 2 0, 

and hence, = gk/[tujJk + J-K" cos(20)/4]. Further, for 
an Ohmic heat bath, the spectral density is given by, 



(59) 



where a is the coupling strength and fi is the cut-off 
frequency of the bath. The renormalization factor K can 
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be obtained from the self-consistent equation, 



k = { 



JKcos(28)/4 
Q + JK cos(2<9)/4 



)«exp 



n + JK cos(26>) /4 



(60) 

Like in the spin-boson model, the parameter K shows 
a crossover behavior below a ~ 1. Two spins become 
parallel to x-direction and the magnetization along x- 
axis sharply becomes 1 when a crosses a critical value a c 
for finite field T. 

Also we focus on another aspect of this wavefunction. 
For 8 — 0, the wave- function |— ) is maximally entangled 
as in a 'cat-state—. In the other limit 8 = ir/4 (in the 
limit of large field T), the ground state |— ) can be factor- 
ized. The amount of entanglement between the two spins 
can be quantitatively calculated from the 'concurrence' 
of the wave function of two spins, which is cos(26>) 24 . 
It is interesting to notice that the concurrence shows a 
crossover at T = by tuning the dissipative coupling 
strength a as shown in Fig. 3. 



VI. SUMMARY AND OUTLOOK 

In this paper we analyze the effect of dissipation on a 
model quantum system which undergoes a QPT. In vari- 
ous limiting cases this model represents different physical 
systems, starting from a TIM in the thermodynamic limit 
(when the number of spins N — > oo) to a generalized spin- 
boson model of two interacting spins (for N — 2). Also 
this model can be viewed as an interacting version of 
a multi-mode Dicke model describing the atom-photon 
interaction. In the context of a large spin our model 
has application to the area of nanomagnets. We have 
developed a self-consistent variational method to study 
the system for any number of spins, and the spin system 
can be described by an effective anisotropic Heisenberg 



model after integrating out the bosonic modes. The ef- 
fective Heisenberg model shows two interesting effects, 
the renormalization of the orginal coupling between the 
z-component of spins similar to the spin-boson model 
and generation of extra terms which couple the x and 
y-components and of spins. In the thermodynamic limit 
(N — > oo) the model represents a TIM in presence of dis- 
sipation and the phase diagram has been obtained within 
a variational-MFA which is in accordance with the sad- 
dle point approximation. It is interesting to note that 
for large N, a new coupling between the x component of 
spins plays an important role and destroys the ferromag- 
netic phase for a smaller value of transverse field. The 
QPT of the model at zero temperature has been studied 
using the Holstein-Primakoff transformation. The exci- 
tation energy vanishes as (r — Tc) 1 / 2 at the quantum 
critical point. 

In the other limit, for N = 2, the model generalizes 
the usual spin-boson model by including a spin-spin in- 
teraction along the z axis. Unlike a TIM, this two-spin 
system shows zero magnetization along the z-axis but dis- 
plays a sharp change in magnetization along the x-axis. 
For small N, strong renormalization of spin-spin coupling 
of the original Hamiltonian becomes important for the 
crossover phenomenon similar to the original spin-boson 
model. From the point of view of quantum information 
theory the ground state of two-spins shows a transition 
from a maximally entangled 'cat-state' to a factorizable 
state where both spins are aligned along the x-axis. 

In conclusion we have studied the effect of dissipation 
on a simple model of an interacting spin system which 
undergoes QPT. Further extension of the model for an- 
alyzing the effect of dissipation in dynamical quenching 
across a quantum critical point that is also relevant for 
the dynamics and thermodynamics of nanomagnets, as 
alluded to in sec. I, is left for future work2£. 
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